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Heretofore neglected decay modes of heavy MSSM Higgs 
bosons into a variety of neutralino pairs may push the LHC 
discovery reach for these crucial elements of an extended Higgs 
sector to nearly the TeV-scale — if sparticle-sector MSSM 
input parameters are favorable. This is well into the so-called 
decoupling region [1], including moderate to low tan/3 values, 
where no known signals exist for said heavy Higgs bosons via 
decays involving solely SM daughter particles, and the lighter 
/i" mimics the lone SM Higgs boson. While the expanse of the 
Higgs to sparticle discovery region is sensitive to a number 
of MSSM parameters, including in particular those for the 
sleptons, its presence is primarily linked to the gaugino inputs 
— in fact, to just one parameter, M2, if gaugino unification 
is invoked. Thus consideration of high vs. low M2 realms 
in the MSSM should be placed on a par with the extensive 
consideration already given to high vs. low tan (5 regimes. 
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Much effort has been devoted to finding Large Hadron 
Collider (LHC) signals for the quintet of Higgs bosons 
expected in the minimal supersymmetric standard model 
(MSSM) — the hghter CP-even h° , the heavier CP-even 
the CP-odd and the charged pair — and 
to delineating the reaches of the various signal modes 
in terms of MSSM input parameters. Key among these 
have naturally been inputs affecting the tree-level Higgs 
sector, Ad A and tan/3, and by far the best-known depic- 
tion of MSSM Higgs boson search projections is the plot 
of the discovery regions (d.r.'s) in the plane formed by 
Ma and tan (3. Most studies of such signatures, including 
all those shown in well-known renditions of this plot, ne- 
glect or exclude the possibility of Higgs boson decays into 
sparticles. The few studies that have investigated such 
sparticle decay modes [2,3] have focused almost exclu- 
sively on the four-lepton plus missing transverse energy 
(i;^?'"'') signature 

pp^ij",A"^X^x!]-4£ + 2x! (1) 

where £ — e,fj, and the Xi's, the lightest supersymmetric 
particles (LSPs), are stable and generate the E'^'^*'. Pre- 
vious studies have also been restrictive in how the x^'s 
decayed, assuming three-body decay modes with virtual 
Z^'s or charged sleptons while neglecting the possibility 
of on-shell two-body decays. 

The MSSM sparticle spectrum is certainly far richer 
in possibilities for Higgs boson decays, many of which 



can also yield the same M + iJ™"**^ signature. And, as 
the masses of the Higgs bosons grow larger, ignoring 
such potential daughter states becomes increasingly ill- 
conceived. The present study is a first-ever serious ex- 
amination of all possible decay chains leading from the 
production of a heavy neutral MSSM Higgs boson, 
or A^, to a 4£ -I- E'rP'''^ final state [4]: 

pp-.H°,A°~. xtXb , X°X° - 4£ + 2x5 (2) 

(a, 6 = 1,2, i,j — 1,2,3,4). In particular, decays in- 
cluding the heavier neutralino (and, to a lesser extent, 
chargino) states beyond X2 are found to dominate the 
signal over significant swaths of the parameter space. 

Naturally, no such signal would exist if the neutralino 
(and chargino) masses were not light enough for the 
Higgs bosons to decay into pairs of them. Yet for the 
heavy ( ^ 500 GeV) and being considered here 
it is in fact un-natural for such decay modes to be ex- 
cluded. To address this question more quantitatively, 
the MSSM input parameters to the neutralino spec- 
trum must be incorporated. These are, at tree-level (to 
which radiative corrections are very minor [5]): tan/3, 
the higgsino mixing parameter /i, and the soft super- 
symmetry (SUSY)-breaking SU{2)^ and U{1)y gaug- 
ino masses, M2 and Mi, respectively. Gaugino unifi- 
cation, which draws support from the projected merg- 
ing of the SU{2)^ and U{1)y coupling strengths at the 
GUT scale in MSSM scenarios, imposes the TeV-scale 
relation Mi ~ ^ta,n'^9wM2 [6]. Fig. 1 shows the inclu- 
sive production rate for Ai + E^^^^ events from LHC i?" 
and production in the plane formed by /i and M2, 
for representative values of Ma and tan (3, and assum- 
ing an integrated luminosity of 300 fb~^. The shape of 
such contours were found to be quite insensitive to val- 
ues of Ad A ^ 400 GeV and tan /3 ?J 5, though rates tended 
to grow with decreasing AIa and increasing tan/3. If M2 
is not too large, dozens to hundreds of signal events are 
expected for all values of though lower \fj,\ values are 
certainly more favorable. Lower values of also en- 
able large event rates for larger values of M2 as shown. 
Also depicted is the percentage of signal events stem- 
ming from , A'^ — > X2X2 decays only. For values of 
l^i] ^ 400 GeV, this generally exceeds 90%. However, for 
lower values of and low to moderate M2 values, decay 
modes featuring the other heavier neutralinos, X3 and/or 
xl, become dominant. The maxima for the signal event 
rate move deeper and deeper into this previously unana- 
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lyzed zone as Ma grows large. Thus the optimal low 
region was overlooked by previous studies limited to only 
the H^,A^ X2X2 decay mode. 

As the Higgs boson masses are raised and heavier neu- 
tralino states become kinematically accessible as decay 
products, it becomes ever more feasible to have slep- 
ton masses situated below the heavier neutralino masses. 
Neutralino decays to on-shell charged sleptons strongly 
enhance branching ratios (BRs) to Icptonic final states, 
with virtually all such decays yielding leptons (which 
tend to be stiffer due to the heavier parent neutralino). 
Thus the heavier ncutralinos may well differ significantly 
from X2 their capacities as conduits for the production 
of charged leptons from heavy Higgs boson production. 
Here though spoiler decay modes via on-shell sneutri- 
nos must not dominate, since these lead to uninteresting 
neutrino-bearing final states. Fortunately this happens 
along relatively narrow strips in parameter space, such 
as when [7] m~ < m~o < m~, a condition though that is 
quite compatible with early LHC sparticle detection. In 
this case the previously-studied X2X2 decay modes yield 
almost no signal events, while at the same parameter 
space point huge event rates from the thusfar neglected 
heavier neutralino channels are quite possible. 

Full event generator-level studies including the simula- 
tion of a typical LHC detector environment (using private 
codes checked against results in the literature) were per- 
formed using HERWIG 6.5 [8] (which obtains its MSSM 
inputs from ISASUSY 7.58 [9]) with the CTEQ6L [10] set 
of parton distribution functions. Higgs bosons were gen- 
erated through the sub-processes qq,gg H^,A°, con- 
servatively ignoring minor additional contributions from 
2^3 subprocesses which would in any event tend to fail 
cuts imposed later. (The same HERWIG processes were 
used to normalize rates shown in FIG. 1.) Four-lepton 
events are first selected according to these criteria: 

0. Events have exactly four leptons, £ = e or /i, consisting 
of two opposite-sign, same-flavor lepton pairs, meeting 
the following criteria: (i) each lepton must have \r]^\ < 2.4 
and > 7, 4 GeV for e, /i, respectively; (ii) each lepton 
must be isolated, set by the demands that there be no 
tracks (of charged particles) with pt > 1.5 GeV in a cone 
of r = 0.3 radians around a specific lepton, and also that 
the energy deposited in the electromagnetic calorimeter 
be less than 3 GeV for 0.05 radians < r < 0.3 radians. 
Selected events are then subjected to the following cuts: 

1. no opposite-charge same-flavor lepton pairs may re- 
construct Mz± 10 GeV. 

2. all leptons must have 20 GeV < £;f < 80 GeV. 

3. events must have 20 GeV < E^'^^ < 130 GeV. 

4. any jets (defined for Et > 20 GeV, I77I < 2.4 and with 
AR = 0.5 [11]) must have < 50 GeV. 

An additional cut on the four-lepton invariant mass was 
also examined; however, the optimal numerical value for 
this cut depends on the masses of the Higgs bosons and 
the LSP, which are assumed to be a priori unknown. If 
these masses are known, even approximately, then this 



extra cut, not employed here, can certainly further re- 
duce the background rates. 

All noteworthy background processes were simulated, 
including those from Standard Model (SM) processes 
Z°Z°^*'> and and MSSM ones involving 7 and/or 

i7 pair production, xx or q/gx pair production (here 
the neutralinos and/or charginos are produced via an s- 
channel off-shell Z^* or ly^*), g and/or g pair produc- 
tion (here leptons are generated via subsequent cascade 
decays featuring charginos and neutralinos, see [12], but 
the jet cut is eff'ective at eliminating such events, assum- 
ing reasonably heavy colored sparticles with masses well 
above those of the colorless EW states), tth°, and tH~ 
+ c.c. (this last process, studied in [13], could arguably 
be lumped together with the signal processes; however, 
the cuts employed here are not designed to select such 
events, with at best only a handful of these surviving). 

Information gathered from FIG. 1 enables selection of 
representative MSSM input parameter sets, aside from 
Ma and tan/3, for the uncolored sparticles (colored 
gluinos and squarks are fixed around the TeV scale). For 
this work, one set is chosen to have 41 + E^^^^ events 
coming almost exclusively from H^, A'^ — > X2X2 decays: 
Set 1: M2 = 180 GeV, Mi = 90 GeV, = -500 GeV, 
m~ = m~ = 250 GeV; 

(soft -^soft ' 

while for a second set channels involving the heavier neu- 
tralinos tend to dominate: 

Set 2: M2 = 200 GeV, Mi = 100 GeV, /x = -200 GeV, 
m~ =150GeV, m~ = 250 GeV. 

Here m~ and m~ ^ are SUSY-breaking slepton mass 

inputs, with m~ = m~ = m~ , to which the physical 

'-■soft 

slepton masses are closely linked [12], and Ar = A^ = 0. 
Signal(S) and background(B) rates for each parameter 
set are studied across the {Ma, tan/3) plane to map out 
potential d.r.'s. The exact criteria used for demarcating 
d.r.'s are presence of at least 10 S events and a 99%- 
confidencc-lcvel upper limit on B events below the 99%- 
confidence-level lower limit on S-l-B events [14]. (As the 
number of S events grows large, this merges with the 
typical S / ^/B ^ 4-5 criterion.) Simulation runs at well 
over a thousand points in the {Ma, tan/3) plane yield 
the results shown in FIG. 2 and FIG. 3 (detailed tabular 
results from representative simulation runs will be pre- 
sented elsewhere) for Set 1 and Set 2, respectively. The 
U + E^'^"" d.r.'s for LHC luminosities of Lint = 100 fb"^ 
and 300 fb~^ are delineated by the thickened contours, 
with the 100 fb~^ contours inscribed within the 300 fb~^ 
ones. Also shown are the region excluded by LEP results 
and the the expected reaches, assuming dnt = 300 fb~^, 
of Higgs boson decay modes into SM daughter particles 
as developed by the ATLAS collaboration [15]. The new 
4:£+E^^^^ d.r.'s clearly provide significant coverage of the 
moderate tan /3, high Ma so-called 'decoupling region' [1] 
inaccessible to the other signals [15-17]. Certainly, the 
Higgs boson to SM daghters' d.r.'s taken from the ATLAS 
results are not obtained using the same choice of MSSM 
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input parameters as are the A° XiX'j ones devel- 
oped in the present work. In fact, the former d.r.'s use 
input choices designed to ehminate, or at least minimize, 
Higgs boson decays into sparticles. This means that, were 
the ATLAS reaches to be generated for the same set of 
neutralino input parameters as the H°,A° XiX^ dis- 
covery regions, the former may well shrink in size, further 
emphasizing the importance of thoroughly studying the 
jjO y^Q ^ 4£ _|_ ijmiss signature. The diminution of the 
expected signatures from SM decay modes of the MSSM 
Higgs bosons was investigated in [2,17]. Assumptions in- 
herent in the ATLAS d.r.'s for the SM decay modes of 
the MSSM Higgs bosons are no less restrictive than the 
choices of MSSM input parameters made to generate the 
two U + E^'^"" d.r.'s in this work. 

Comparing d.r.'s for Set 1 and Set 2, startling differ- 
eces are immediately apparent (these will be described in 
more detail elsewhere [18]). Of perhaps foremost impor- 
tance among these is Set 2's far greater stretch to high 
Ma values. This is due to the heavier neutralino decay 
modes first described herein: in FIG. 3 corresponding to 
Set 2, inclusion of only X2X2 decay modes, as in pre- 
vious studies, yields a mere smidgen of a d.r. aroimd ( 
Ma, tan/3 ) = ( 350-400 GcV, 3-4.5 ) rather than the ex- 
pansive region shown. All of the large portion of the d.r. 
above tan/3 = 5 would be lost. Note how there is some 
d.r. reach for almost all possible values of tan /3, demon- 
strating the relative merit of the M2 and /z parameters in 
demarcating heavy Higgs bosons' d.r.'s. Previously, par- 
ticular attention to the high tan (3 regime of the MSSM 
has been afforded in this journal [19]. It should however 
be noted that the parameter choices considered herein 
are no more limited in scope than studies devoted exclu- 
sively to that subset of model parameter space. Indeed, 
artificially foisting sparticle masses into the stratosphere 
to avoid their presence in Higgs boson decays, as is the 
norm in most MSSM Higgs boson studies (which then 
resemble two-Higgs doublet model studies) is arguably 
even more restrictive. For all their differences, roughly 
representing different extremes in the signal composition, 
both Set 1 and Set 2 offer new coverage for the crucial 
moderate to high Ma, middle tan/3 zone. 

In conclusion, restricting the search for heavy MSSM 
Higgs boson LHC signals to SM decay modes plus per- 
haps the single neutralino pair X2X2 misses a significant 
portion of the potential d.r. accessible with the inclusion 
of all sparticle decays. A complete analysis featuring all 
sparticle channels yielding the M + E!^^^^ signature is 
quite tractable, and, for the heavy iJ" and A'^ present 
in the so-called decoupling region, (i) this is the only 
known signal mode; and (ii) most of the discovery re- 
gion mapped out relies on the heavier neutralinos, which 
also more easily accommodate rate-enhancing underlying 
sleptonic states, never included in such analyses before. 
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FIG. 1. Number of pp H", A" AIN events, where £ — or /i* and A'' represents invisible final state particles, for 
300fb~\ with tan/3 = 10 and Ma = 500 GeV. Also shown is the percentage from H°,A° X2X2 (in color / black & white): 
> 90% (red cross-hatched / grey cross-hatched), 50% - 90% (red/dark shaded), 10% - 50% (blue/lighter shading), < 10% 
(white/no shading). Optimized slepton masses [13] (with stau inputs raised 100 GeV) are used, rrit = 175 GeV, mt — 4.25 GeV, 
m~= 1 TeV, m~ = 800 GeV, Ar — Ap = 0. The blackened-out areas are excluded by LEP. 
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FIG. 2. Discovery region in (MA,tan/3) plane for MSSM Parameter Set 1 and dm = 100 fb"^ & 300 fb"^ for MSSM Higgs 
bosons' 4£ signals from their XiX^jXmXn decays (here H^jA^ decays to X2X2 totally dominate), shown together with regions 
for other MSSM Higgs boson signatures from decays to SM particles based upon LEP results and ATLAS simulations [15] 
(which assume dnt ~ 300 fb~^) where labels represent: 1. — > Z'^ Z'^* — > 41eptons; 2. t — > hH^ , t^v -\- c.c; 3. 

tih°, h° ~* bb; 4. h° 77 and W^h°/tth°, /i° ^ 77; 5. bbH°, bbA° with 7^7^° bb; 6. H+ ^ tb + c.c; 7. H°/A° fi+ n' ■ 
8. H°/A° T+r-; 9. gb tH+, H+ t+v + c.c; 10. H° h°h° bb'yy, 11. A° Z°h^' ~* l+Tbl; 12. H° /A'' ti. 
Note that SM discovery regions are not for the same input parameters: they presume Higgs bosons cannot decay into sparticles, 
so more accurate estimates may well be smaller. 

ft 
as 




800 
(GeV) 

FIG. 3. Discovery region in (Ma, tan/3) plane for for MSSM Parameter Set 2 and dnt = 100 fb"^ & 300 fb"^ for MSSM 
Higgs bosons' 4£ signals from their XiXji XmXn. decays (here Higgs boson decays to higher-mass neutralinos typically dominate), 
shown together with regions for MSSM Higgs boson signatures from decays to SM particles as in FIG. 2. 
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